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Abstract 

The problem of finding of the ferromagnetic and antiferromagnetic "symmetry broken" 
solutions of the correlated lattice fermion models beyond the mean-held approximation has 
been investigated. The calculation of the quasiparticle excitation spectra with damping 
for the single- and multi-orbital Hubbard model has been performed in the framework of 
the equation- of-motion method for two-time temperature Green's Functions within a non- 
perturbative approach. A unified scheme for the construction of Generalized Mean Fields 
(elastic scattering corrections) and self-energy (inelastic scattering) in terms of the Dyson 
equation has been generalized in order to include the presence of the " source fields" . The 
damping of quasiparticles, which reflects the interaction of the single-particle and collective 
degrees of freedom has been calculated. The "symmetry broken" dynamical solutions of 
the Hubbard model, which correspond to various types of itinerant antiferromagnetism has 
been discussed. This approach complements previous studies and clarifies the nature of 
the concepts of itinerant antiferromagnetism and "spin-aligning field" of correlated lattice 
fermions. 
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1 Introduction 



The problem of the adequate description of strongly correlated lattice fermions has been studied 
intensively during the last decade, especially in the context of Heavy Fermions and High-Tc 
superconductivity [U - |||. The behaviour and the true nature of the electronic states and 
their quasiparticle dynamics is of central importance to the understanding of the magnetism in 
metals and the Mott-Hubbard metal-insulator transition in oxides, the heavy fermions in rare- 
earths compounds and the high-temperature superconductivity (HTSC) in cuprates. Recently 
there has been considerable interest in identifying the microscopic origin of these states Q. 
Antiferromagnetic correlations may play an important role in the possible scenario of normal 
and superconducting behavior of these compounds. Some of the experimental and theoretical 
results show that antiferromagnetic spin fluctuations are really involved in the problem. This 
idea has stirred a great deal of discussion in recent times ||. An appealing but phenomenological 
picture of HTSC, known as the nearly antiferromagnetic Fermi liquids (NAFL) approach, has 
been developed to explain many anomalous properties of cuprates M. This approach predicts the 
detailed phase diagram for cuprates || and present arguments which suggest that the physical 
origin of the pseudogap found in quasiparticle spectrum below the critical temperature is the 
formation of a precursor to a spin-density-wave-state. While the NAFL's scenario is appealing, 
it has aparently not yet been derived from fully microscopic consideration. The problem of 
the role of antiferromagnetic spin fluctuations for HTSC has recently been the subject of many 
papers( for a recent review see e.g. Ref. jjj). These investigations call for a better understanding 
of the nature of solutions (especially magnetic) to the Hubbard and related correlated models [|] 
- 1 11]. The microscopic theory of the itinerant ferromagnetism and antiferromagnetism [12], ]l3| 
of strongly correlated fermions on a lattice at finite temperatures is one of the important issues 
of recent efforts in the field [14] - 17]. In some papers the spin-density- wave (SDW) spectrum 
was only used without careful and complete analysis of the quasiparticle spectra of correlated 
lattice fermions. The aim of this paper is to investigate the intrinsic nature of the "symmetry 
broken" ( ferro- and antiferromagnetic) solutions of the Hubbard model at finite temperatures 
from the many-body point of view. In the previous papers we set up the formalism and derived 
the equations for the quasiparticle spectra with damping within the single- and multi-orbital 
Hubbard model for the uniform paramagnetic case. In this paper we apply the formalism to 
consider the ferromagnetic and antiferromagnetic solutions. It is the purpose of this paper to 
explore more fully the notion of Generalized Mean Fields (GMF) []l0[ which may arise in the 
system of correlated lattice fermions to justify and understand the "nature" of the local staggered 
mean-fields which fix the antiferromagnetic ordering. The present work brings together the 
formulation of the itinerant antiferromagnetism of various papers. For this aim we rederive the 
SDW spectra by the Irreducible Green's Functions (IGF) method taking into account the 
damping of quasiparticles. This alternative derivation has a close resemblance to that of the 
BCS theory of superconductivity for transition metals [19], [20] using the Nambu representation 
(c.f. |2^]).This aspect of the theory is connected with the concept of broken symmetry, which is 
discussed in detail for the present case. The advantage of the Green's function method is the 
relative ease with which temperature effects may be calculated. 
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2 Itinerant Antiferromagnetism 



The antiferromagnetic state is characterized by a spatially changing component of magnetisation 
which varies in such a way that the net magnetisation of the system is zero. The concept of 
antiferromagnetism of localized spins which is based on the Heisenberg model and the two- 
sublattice Neel ground state is relatively well founded contrary to the antiferromagnetism of 
delocalized or itinerant electrons. The itinerant-electron picture is the alternative conceptual 
picture for magnetism 22]. 



We now sketch the main ideas of the concept of itinerant antiferromagnetism. The simplified 
band model of an antiferromagnet has been formulated by Slater [p3] within the single-particle 
Hartree-Fock (H-F) approximation. In this approach he used the "exchange repulsion" to keep 
electrons with parallel spins away from each other and to lower the Coulomb interaction energy. 
Some authors consider it as a prototype of the Hubbard model. However the exchange repulsion 
was taken proportional to the number of electrons with the same spins only and the energy 
gap between two subbands was proportional to the difference of electrons with up and down 
spins. In the antiferromagnetic many-body problem there is an additional "symmetry broken" 
aspect. For an antiferromagnet, contrary to ferromagnet, the one-electron H-F potential can 
violate the translational crystal symmetry. The period of the antiferromagnetic spin structure 
L is greater than the lattice constant a. To introduce the two-sublattice picture for itinerant 
model one should assume that L = 2a and that the spins of outer electrons on neighboring 
atoms are antiparallel to each other. In other words, the alternating (H-F) potential Vi a = 
— av ex.p(iQRi) where Q = (vr/2, ir/2, ir/2) corresponds to a two-sublattice AFM structure. To 
justify an antiferromagnetic ordering with alternating up and down spin structure we must admit 
that in effect two different charge distributions will arise concentrated on atoms of sublattices 
A and B. This picture accounts well for quasi-localized magnetic behavior. 
The earlier theories of itinerant antiferromagnetism were proposed by des Cloizeaux ]24| and 
especially Overhauser (2[| (in the context of the investigation of the ground state of nuclear 



matter). Then Overhauser [26] has applied this approach for the explanation of the anomalous 
properties of dilute Cu — Mn alloys, has suggested an antiferromagnetic mechanism that requires 
neither two-body interactions between paramagnetic solute spins, nor a sublattice structure 
(c.f. [0]). Such a mechanism may be recognized by considering a new type excited state of 
the conduction electron gas. He invented the static SDW which allow the total charge density 



of the gas to remain spatially uniform. Overhauser pi] - [29] suggested that the H-F ground 



state of a three dimensional electron gas is not necessarily a Slater determinant of plane waves. 
Alternative sets of one-particle states can lead to a lower ground-state energy. Among these 
alternatives to the plane-wave state are the SDW and CDW ground states for which the one- 
electron Hamiltonians have the form 

H = (p 2 /2m) — G(a x cos Qz + a y sinQz) (1) 

( spiral SDW; Q = 2k F z ) 
and 

H = (//2m) - 2G cos(Qr) (2) 

( CDW; Q = 2k F z ) 

The periodic potentials in the above expressions lead to a corresponding variation in the elec- 
tronic spin and charge densities, accompanied by a compensating variation of the background. 
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The effect of Coulomb interaction on the magnetic properties of the electron gas in Overhauser's 
approach renders the paramagnetic plane-wave state of the free-electron-gas model unstable 
within the H-F approximation. The long-range components of the Coulomb interaction are 



most important in creating this instability []2S|]. It was demonstrated 28] that a nonuniform 
static SDW is lower in energy than the uniform (paramagnetic state) in the Coulomb gas within 
the H-F approximation for certain electron density. 

The H-F is the simplest approximation but neglects the important dynamical part. To include 
the dynamics one should take into consideration the correlation effects. The role of correlation 
corrections which tend to suppress the SDW state as well as the role of shielding and screen- 
ing were not fully clarified |3C]. Overhauser remarked that SDW ground states do not occur 



for (^-function interactions, whatever their strength. This question was investigated further in 



Ref. [31|. An instability of the paramagnetic Hartree-Fock state against a state with different 
orbitals for different spins was interpreted as a magnetic phase transition. 

It is important to note that in the Slater's and des Cloiseaux's models an electron moving in a 
crystal does not change its spin. In these models the main processes are related with the pairing 
of electrons having the same spins, one from each of the two sublattices. In the Overhauser's 
approach to itinerant antiferromagnetism the combination of the electronic states with different 
spins (with pairing of the opposite spins) is used to describe the SDW state with period Q. The 
first approach is obviously valid only in the simple commensurate two-sublattice case and the 
latter is applicable to the more general case of an incommensurate spiral spin state. The general 
SDW state has the form 

= X P a cos(6> p /2) + Xp+Q-v sin(0 p /2) (3) 

The average spin for helical or spiral spin arrangement changes its direction in the (x-y) plane. 
For the spiral SDW states a spatial variation of magnetisation corresponds to Q = (^)(1, 1)- 
The antiferromagnetic phase of chromium [32], [33] and its alloys has been satisfactorily ex- 



plained in terms of the SDW within a two-band model |3j] . It is essential to note that chromium 
becomes antiferromagnetic in a unique manner. The antiferromagnetism is established in a more 
subtle way from the spins of the itinerant electrons than the magnetism of collective band elec- 
trons in metals like iron and nickel. The essential feature of chromium which makes possible 
the formation of the SDW is the existence of "nested" portions of the Fermi surface (3^]. The 
formation of bound electron-hole pairs takes place between particles of opposite spins; the con- 
densed state exhibits the SDW. 

The recent attempt to describe an antiferromagnetic insulator at T = using a one-electron 
approach was made in Ref. [j35|. To do this, the authors proposed to overcome the inadequacies 
of standard local-spin-density theory by adding a spin-dependent magnetic pseudopotential to 
Kohn-Sham equations. 



For the Hubbard model [36] the qualitative phase diagram was calculated by Penn [37]. Un- 
fortunately, although his work gives a clear physical picture, it does not emphasize the lattice 
character of the tight-binding or Wannier fermions as well as the essence of the anomalous spin- 
flip averages. The Hubbard model is a simplified but workable model for the correlated lattice 
fermions and the applicability of the SDW Overhauser concept to highly correlated tight binding 



electrons on a lattice deserves a careful analysis within this model. In earlier papers [38] - [41] the 
single- and multi-orbital Hubbard model has been investigated with respect to antiferromagnetic 
solutions in the mean-field approximation mainly. 
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3 Hubbard Model 



The Hubbard model has been widely recognised as a workable model for a study of the correlated 
itinerant electron systems. For the sake of completeness we shall discuss the single-orbital and 
multi-orbital cases separately. 

3.1 MULTI-ORBITAL HUBBARD MODEL 

To demonstrate the advantage of our approach we shall consider the quasiparticle spectrum of 
the lattice fermions for degenerate band model. Let us start with the second quantized form of 
the Hamiltonian taking the set of the Wannier functions [4>\{r — Ri)]. Here A is the band index 
(A= 1,2,. ..5). 

H= X t ij a icr a j^ + 2 X X <iatJP\W\mi,nS>a^afp a ,a my<r >a n s ir (4) 

ijfiua ij,mn a/3-ySaa' 

For a degenerate d-band the second quantized form of the total Hamiltonian in the Wannier- 
function representation reduces to the following model Hamiltonian 

H = H 1 +H 2 + H 3 (5) 

The kinetic energy operator is given by 

ffi=EEfiv (6) 

ij ^va 

The term H2 describes one-centre Coulomb interactions 

H 2 = - ^2 U ^ n i^ n i^-o- + ^ X X V ^ n i^ n iy<r' ( l ~ <W) ~ ( 7 ) 

ifii/cr ifii/a 

9 X ^ a i/icr a *M-o- a jt-o- a ^o"(l — Spa,) 
i /ii/CT 

In addition to the intrasite intraorbital interaction which is the only interaction present in 
the single-orbital Hubbard model, this term contains three more kinds of interorbital interac- 
tions. 

The last term H3 describes the direct intersite exchange interaction 

H 3 = - g X X J ^ a t* a in-<T' a tn*i a W° ( 8 ) 

ij^i air' 

The definition of various integrals in H is obvious. It is reasonable to assume that: 

U^ = U; V^ = V; 1^ = 1; Jff = Jy. (9) 



This Hamiltonian differ slightly from the analogous Hamiltonian of Ref. [4C] where the only 
intrasite interaction terms of the second-quantized Hamiltonian of the d-band were taken into 
consideration. 
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3.2 SINGLE-ORBITAL HUBBARD MODEL 



The model Hamiltonian which is usually referred to as Hubbard Hamiltonian [36] 

H = ^2t ij af a a ja + U/2^2n ia ni- a (10) 

ijtr i(r 

includes the intraatomic Coulomb repulsion U and the one-electron hopping energy tij. The 
electron correlation forces electrons to localize in the atomic orbitals, which are modelled here 
by the complete and orthogonal set of the Wannier wave functions [4>{r — Rj)]. (The Wannier 
representation, which is a unitary transformation of the Bloch representation is an important 
background of the Hubbard model. It is well known that in one-dimension the Wannier functions 
decrease exponentially but less is known about two- and three-dimensional cases.) On the other 
hand, the kinetic energy is reduced when electrons are delocalized. The main difficulty of the 
right solution of the Hubbard model is the necessity of taking into account both these effects 
simultaneously. Thus, the Hamiltonian (10) is specified by two parameters: U and effective 
electron bandwidth 

A = (N- 1 Y,\U 3 \ 2 ) 1/2 - 

The important third "player" is the Pauli principle, which has a long-range character, contrary 
to the local Coulomb repulsion and nearest-neighbour hopping. 
The band energy of Bloch electrons e(k) is defined as follows 

t {j = N~ l exp[ik(Ri - Rj)}, 

k 

where ./V is the number of the lattice sites. It is convenient to count the energy from the 
center of gravity of the band, i.e. ta = J2k e (^) = 0- The effective electron bandwidth A and 
Coulomb intrasite integral U define completely the different regimes in 3 dimension depending 
on parameter 7 = A/U. It is usually a rather difficult task to find interpolating solution for 
the dynamical properties of the Hubbard model. We evidently have to to improve the early 
Hubbard's theory taking account of the variety of possible regimes for the model depending on 
electronic density, temperature and values of 7. It was the purpose of the papers ||, |l(| to 
find the electronic quasiparticle spectra in a wide range of the temperature and the parameters 
of the model and to account explicitly for the contribution of damping of the electronic states 
when calculating the various characteristics of the model. In the past years many theoretical 
papers have been published , in which the approximative dynamical solution of the models (5) 
and (10) has been investigated by means of various advanced methods of many-body theory. 
Despite the considerable contributions to the development of the many-body theory and to 
our better understanding of the physics of the correlated electron systems, the fully consistent 
dynamical analytical solution of the Hubbard model is still lacking. To solve this problem 
with a reasonably accuracy and to describe correctly an interpolating solution one need a more 
sophisticated approach than the usual procedures which have been developed for description of 
the interacting electron-gas- problem. 

4 Irreducible Green's Functions Method 

Recent theoretical investigations of strongly correlated electron systems have brought forward 
a significant variety of approaches. To describe from first principles of the condensed matter 
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theory and statistical mechanics the physical properties of strongly correlated systems we need 
to develop a systematic theory of quasiparticle spectra. 

In this paper we will use the approach which allows one to describe completely the quasi-particle 
spectra with damping in a very natural way. This approach has been suggested to be essential for 
various many-body systems and we believe that it bears the real physics of strongly correlated 
electron systems [flOfl , [18]. The essence of our consideration of the dynamical properties of many- 



body system with strong interaction is related closely with the field theoretical approach and 
use the advantage of the Green's functions language and the Dyson equation. It is possible to 
say that our method tend to emphasize the fundamental and central role of the Dyson equation 
for the single-particle dynamics of the many-body systems at finite temperatures. 
In this Section, we will discuss briefly this novel nonperturbative approach for the description 
of the many-body dynamics of strongly correlated systems. A number of other approaches has 
been proposed and our approach is in many respect additional and incorporates the logic of 
development of the many-body techniques. The considerable progress in studying the spectra of 
elementary excitations and thermodynamic properties of many-body systems has been for most 
part due to the development of the temperature dependent Green's Functions methods. We have 
developed a helpful reformulation of the two-time GFs method which is especially adjusted [|| 
for the correlated fermion systems on a lattice. The very important concept of the whole method 
are the Generalized Mean Fields. These GMFs have a complicated structure for the strongly 
correlated case and do not reduce to the functional of the mean densities of the electrons, when 
we calculate excitations spectra at finite temperatures. To clarify the foregoing, let us consider 
the retarded GF of the form 

G r =« A(t), B{t') »= -i0(t - < [A(t)B(t% >, V = ±1. (11) 

As an introduction of the concept of IGFs let us describe the main ideas of this approach in a 
symbolic form. To calculate the retarded GF G(t — t') let us write down the equation of motion 
for it: 

u G{u) =< [A,A + ] V > + « [A, H]^ | A + » w . (12) 

The essence of the method is as follows JL8]. It is based on the notion of the "IRREDUCIBLE" 
parts of GFs (or the irreducible parts of the operators, out of which the GF is constructed) in 
term of which it is possible, without recourse to a truncation of the hierarchy of equations for the 
GFs, to write down the exact Dyson equation and to obtain an exact analytical representation 
for the self-energy operator. By definition we introduce the irreducible part (ir) of the GF 

ir « [A,H]-\A + »=« [A,H]- - zA\A + » . (13) 

The unknown constant z is defined by the condition (or constraint) 

< [[A, H]-, A + ] ri >= (14) 

From the condition (14) one can find: 

_ < [[A,H]„,A+] TI > _ Mi 



<[A,A+] V > M 



(15) 



Here Mq and Mi are the zeroth and first order moments of the spectral density. Therefore, 
irreducible GF are defined so that they cannot be reduced to the lower-order ones by any kind 
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of decoupling. It is worthy to note that the irreducible correlation functions are well known in 
statistical mechanics. In the diagrammatic approach the irreducible vertices are defined as the 
graphs that do not contain inner parts connected by the G°-line. With the aid of the definition 
(13) these concepts are translated into the language of retarded and advanced GFs. This pro- 
cedure extract all relevant (for the problem under consideration) mean field contributions and 
puts them into the generalized mean- field GF, which here are defined as 

G o = <[A,A+] > ^ 
(uj-z) 

To calculate the IGF ir « [A, H]_(t), A + (t') » in (12), we have to write the equation of 
motion after differentiation with respect to the second time variable t' . The condition (14) 
removes the inhomogeneous term from this equation and is a very crucial point of the whole 
approach. If one introduces an irreducible part for the right-hand side operator as discussed 
above for the "left" operator, the equation of motion (12) can be exactly rewritten in the 
following form 

G = G° + G°PG°. (17) 

The scattering operator P is given by 

P=(M )~ 1 ir « [A,H]_\[A + ,H}_ » ir (M )-\ (18) 

The structure of the equation (17) enables us to determine the self-energy operator M, in 
complete analogy with the diagram technique 

P = M + MG°P. (19) 

From the definition (19) it follows that the self-energy operator M is defined as a proper (in 
diagrammatic language "connected") part of the scattering operator M = (P) p - As a result, we 
obtain the exact Dyson equation for the thermodynamic two-time Green's Functions: 

G = G° + G°MG, (20) 

which has a well known formal solution of the form 

G = [(G )- 1 - M}- 1 (21) 

Thus, by introducing irreducible parts of GF (or the irreducible parts of the operators, out 
of which the GF is constructed) the equation of motion (12) for the GF can be exactly (but 
using constraint (14)) transformed into Dyson equation for the two-time thermal GF. This is 
very remarkable result, which deserves underlining, because of the traditional form of the GF 
method did not include this point. The projection operator technique has essentially the same 
philosophy, but with using the constraint (14) in our approach we emphasize the fundamental 
and central role of the Dyson equation for the calculation of the single-particle properties of 
the many-body systems. It is important to note, that for the retarded and advanced GFs the 
notion of the proper part is symbolic in nature fl8fl . However, because of the identical form of 
the equations for the GFs for all three types (advanced, retarded and causal), we can convert 
in each stage of calculations to causal GFs and, thereby, confirm the substantiated nature of 
definition (19)! We therefore should speak of an analogue of the Dyson equation. Hereafter we 
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will drop this stipulation, since it will not cause any misunderstanding. It should be emphasized 
that the scheme presented above give just an general idea of the IGF method. The specific 
method of introducing IGFs depends on the form of operator A, the type of the Hamiltonian 
and the conditions of the problem. The general philosophy of the IGF method lies in the sepa- 
ration and identification of elastic scattering effects and inelastic ones. This last point is quite 
often underestimated and both effects are mixed. However, as far as the right definition of 
quasiparticle damping is concerned, the separation of elastic and inelastic scattering processes 
is believed to be crucially important for the many-body systems with complicated spectra and 
strong interaction. Recently it was emphasized especially that the anomalous damping of elec- 
trons (or holes) distinguishes cuprate superconductors from ordinary metals. From a technical 
point of view the elastic (GMF) renormalizations can exhibit a quite non-trivial structure. To 
obtain this structure correctly, one must construct the full GF from the complete algebra of the 
relevant operators and develop a special projection procedure for higher-order GF in accordance 
with a given algebra. It is necessary to emphasize that that there is an intimate connection 
between the adequate introduction of mean fields and internal symmetries of the Hamiltonian. 



5 Symmetry Broken Solutions 

In many-body interacting systems, the symmetry is important in classifying of the different 
phases and in understanding of the phase transitions between them [42] - |48j. According to 



Bogolubov p2| ( c.f. |47j ) in each condensed phase, in addition to the normal process, there is an 
anomalous process (or processes) which can take place because of the long-range internal field, 
with a corresponding propagator. The anomalous propagators for interacting many-fermion sys- 
tem corresponding to the ferromagnetic (FM) and antiferromagnetic (AFM) long-range ordering 
are given by 

FM : Gf m ~« a ka , a+_ a » (22) 
AFM : G afm ~<< a k+Qa -al +QI(j , » 

In the SDW case, a particle picks up momentum Q — Q' from scattering against the periodic 
structure of the spiral (nonuniform) internal field, and has its spin changed from a to a 1 by the 
spin-aligning character of the internal field. The Long-Range-Order (LRO) parameters are: 

FM :m = l/Nj2< 4* a k-* > (23) 

kc 



AFM :M Q = J2< 4* a k+Q-* > 

ka 

It is important to note that the long-range order parameters are functions of the internal field, 
which is itself a function of the order parameter. There is a more mathematical way of for- 
mulating this assertion. According to the paper [14|, the notion "symmetry breaking" means 



that the state fails to have the symmetry that the Hamiltonian has. True broken symmetry 
can arise only if there are infinitesimal "source fields" present. Indeed, for the rotationally and 
translationally invariant Hamiltonian the suitable source terms should be added: 

FM : efi B H x 4a a k-a (24) 

ka 

AFM : EH B H^2,4a a k+Q-cT 

kQ 
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where e — > is to be taken at the end of calculations. 

Broken symmetry solutions of the Overhauser type (3) imply that the vector Q is a measure of 
the inhomogeneity or breaking of translational symmetry. It is interesting to note the remark of 
paper |iq| (c.f . p8| ) about antiferromagnetism, for which "a staggered magnetic field plays the 
role of symmetry-breaking field. No mechanism can generate a real staggered magnetic field in 
an antiferromagnetic material" . The Hubbard model is a very interesting tool for the analyzing 
this concept [41| 



Penn [37] shown that antiferromagnetic state and more complicated states (e.g. ferrimagnetic) 
can be made eigenfunctions of the self-consistent field equations within an "extended" mean- 
field approach, assuming that the "anomalous" averages < af a ai- a > determine the behavior 
of the system on the same footing as the "normal" density of quasiparticles < af a ai a >. It 
is clear, however, that these "spin-flip" terms break the rotational symmetry of the Hubbard 
Hamiltonian. For the single- band Hubbard Hamiltonian the averaging < af_ a ai t(T >= because 
of the rotational symmetry of the Hubbard model. The inclusion of the "anomalous" averages 
lead to the unresricted H-F approximation. The rigorous definition of the unrestricted Hartree- 



Fock approximation (UHFA) has been done recently in Ref. [14|. This approximation has been 
applied also for the single-band Hubbard model (10) for the calculation of the density of states. 
The following definition of UHFA has been used: 

ni- a a ir7 =< rii^cr > a ir7 - < al_ a a ia > ai- a (25) 

Thus, in addition to the standard H-F term, the new, the so-called "spin-flip" terms, are retained. 
This example clearly show that the nature of the mean-fields follows from the essentials of the 
problem and should be defined in a proper way. So, one needs a properly defined effective 
Hamiltonian H c r . We shall analyze below in detail the proper definition of the irreducible GFs 
which include the "spin-flip" terms. For the single-orbital Hubbard model this definition should 
be modified in the following way: 

iv + I + + I + 

<< a>kJ r p a (lpj r q_ a O,q— o^kfj ^ > ^ > ui = < ^ < ^ a 'k+p<jQ'p-i r q—<jQ j q—o\Q'kcF — 

<5p,0 < n q ^ a > Gka~ < dk+paap+q-o- ><< a q-a\ a ka »u (26) 

From this definition it follows that such way of introduction of the IGF broadens the initial 
algebra of the operators and the initial set of the GFs. This means that "actual" algebra of the 
operators must include the spin-flip terms at the beginning, namely: (aj CT , af a , m a , a^ T aj_ (J ). 
The corresponding initial GF will have the form 



« a irT \aJ_ c » 



« ai- a \aj a » « cn- a \aj_ a » 



With this definition we introduce the so-called anomalous (off-diagonal) GFs which fix the 
relevant vacuum and select the proper symmetry broken solutions. In fact, this approximation 



has been investigated earlier by Kishore and Joshi 5Cj. They clearly pointed out that they 
assumed that the system is magnetized in the x direction instead of the conventional z axis. 
The detailed investigation and classification of the magnetic and non-magnetic symmetry broken 
solutions of the three-band extended Hubbard model for Cu02 planes of high- T c superconductors 



was made in Ref. [51] within the mean-field approximation. 
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6 Dynamical Properties 



In many-body interacting systems the quasiparticle dynamics can be quite non-trivial. Here the 
problem of the adequate description of the many-body dynamics of the multi-orbital Hubbard 
model will be discussed in the framework of the equation-of-motion approach for two-time ther- 
modynamic Green's Functions. Our main motivation was the intention to formulate a consistent 
theory of dynamical properties of the Hubbard model taking into account the symmetry broken 
(magnetic) solutions. 

This formulation gives us an opportunity to emphasize some important issues about the relevant 
dynamical solutions of the strongly correlated models of fermions on a lattice and to formulate 
in a more sharp form the ideas of the method of the Irreducible Green's Functions (IGF) [18]. 
This IGF method allows one to describe the quasiparticle spectra with damping of the strongly 
correlated electron systems in a very general and natural way and to construct the relevant dy- 
namical solution in a self-consistent way on the level of the Dyson equation without decoupling 
the chain of the equation of motion for the GFs. 

The interplay and the competition of the kinetic energy and potential energy affects substan- 
tially the electronic spectrum. The renormalized electron energies are temperature dependent 
and the electronic states have finite life times. These effects are most suitably accounted for by 
the Green's functions method. We shall use the (IGF) method of Section 4. To give a more 
instructive discussion let us consider the single- particle GF of lattice fermions, which is defined 
as 

GZSr,t- =« a itl( r(t),af a ,(t') »= -iO{t - t') < [a iflcr (t), a+ a , (*')]+ >= 

/+oo 
duexp(-iut)G%,(ij;u) (27) 
-oo 

Actually, this GF is a matrix (10x10) in the joint tensor product vector space of spin and orbital 
momentum. The diagonal elements of this matrix GF are normal propagators, while the off- 
diagonal elements are anomalous. The equation of motion for the Fourier transform of the GF 
has the form 

]T A» a (im)GZ> (mj; u) = 5^5^, + £[5^) « a miMT n maa \af va , » (28) 

ma ma 

+B% a (im) « a mfl(T n ma ^ a \a+ u(7 , » + 
B% a (im)(« a iflcr n mfia \a+ ual » + « a lfl - a a+ lfl _ a a miI(7 \a+ U(Tl >>)] 

Here we have introduced the notations 



A» a (im) = u5 mi 8 m - C; B^{im) = (V - I)5 im (l - 8^); (29) 
B% a = [US^a + V(l - <V*Rm; B% a (im) = J im (l - 5j m )<V* 

Let us introduce, by definition, an "irreducible" GF in the following way 

( << CLi/3<T a ma(T 1 a maai\o^ UIJ i >>) =<< a i/3cr a mocai a maai \Oj v(r i >> (30) 
— < n maai > m i << a>if) a \ a j va i < ( ^ificfO'maa\ X< dmaai V^jvo' 

According to (14), the following constraint should be valid 

< [(aif3*n maai )( ir \a+ ual } + >= (31) 



10 



Substituting (30) in (28) we obtain the following equation of motion in the matrix (in spin space) 
form 

Y,F^(im)G^(mj;u>) = 1 + ^H^C'Vi) + (32) 

ma ma 

L^(im)D^(mj)+L^(im)(Rr(im,j) + ^(im,j))} 



where 



L^fzm) - ( B H im) TV L m (im) — ( ° ° 

o o J' L2W- Ur(im) srw 

^ (?mj " v B% a (im) 



(33) 



and 



Eg(im) = A^(im) - B^(im) < a m ^a+ a ^ > - (34) 

J2(B^(im) < n maT > 5^ - 
P 

-B% P (im) < n ial > 5^) - 

B% a (im)(< a iaT a+ aT > + < a ial a^ al >)-^B%*(ml) < n M >; 

l 

E^2 = ~ B T < a m rf a mal > ~J2 B 3 a < a ta[ a l^ > $im (35) 

I 

and similar expressions for #21 and £22 with reversed spin indices. The higher-order GF have 
the form 

_ ( ( (ir) << a m/lT n maT |at T >>) (("*)<< a mjuT n m/xT |at { »)\ 
1 \ y ( (ir) << a m ^n rnal \a^ >>) ( (ir) << a m ^n mal \a^ ul ») J 
= ( ( (ir) « a ia] n mar \a^ >>) ( (ir ^« a ia] n ma] |a+ i »)\ 
2 \ ((* r )<< a iai n ma Ja+ T >>) (( ir )<< a iai n mai [0+^ >>) / 



(36) 



and R has the following structure 
R = 



( { (ir) « am^a+ a ^a ma ^\a^ ») ((* r )« a m ^a^a ma] -\a+ H »)\ 
V ( (tr) « am^a+ al a mal \a+ u] ») ((* r )« am M T a mai a ^ll a n^ »)/ 

To calculate the higher-order GF D±, D2, R\ and i?2, we will differentiate the r.h.s. of it with 
respect to the second-time variable (f). Combining both (the first- and second-time differen- 
tiated) equations of motion we get the " exact" ( no approximation have been made till now) 
" scattering" equation 

CT{ir,u) = G^{ij-u) + ]T G^(im;u)P^(mn;uj)G^(nj;u) (37) 

mna/3 
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Here we have introduced the generalized mean-field (GMF) GF Go according to the following 
definition 

J2F^(im)Gr(mj;co) = 5^ (38) 

ma 

The scattering operator P has the form 

P^imn-u) = ( %f> n > U \ %\ mn ' U \) ; (39) 
Let us write down explicitly the first matrix element 

P°f(mn;u) = E[5f°N( W « a m ^n ma ^n m »^)Bf(nj) + (40) 

Br(im)(W« a m ^n ma ^n n(n »^)B^(nj) + 
CW( (!r) « a mil} n mai \a^n m »^)B? v (nj) + 



B^ a (im)(^« a m ^n mai \a+ n m »^)B^(nj)) 



Here we presented for brevity the explicit expression for a part of Hamiltonian (5) only without 
the last term. Using (17) - (19) we find the Dyson equation in the Wannier basis 

G^(ij;uj) = G^(ij;u)+ £ G^(im;u)M a f(mn;uj)G^(nj;u) (41) 

mna/3 

The equation (41) is the central result of the present treatment. 

7 Quasiparticle Formulation 

Let us first consider how to describe our system in terms of quasiparticles. For a translationally 
invariant system, to describe the low-lying excitations in terms of quasiparticles one has to make 
a Fourier transformation 

CT{ir, oj) = N- 1 exp[ik(Ri - Rj)]G tiV (k; u) (42) 

k 

M^(ij; uj) = N- 1 exp[ik(Ri - Rj)]M^{k- u) 
k 

t w = N -i J- exp[ik(Ri - Rj)]e^k) 
k 

The Dyson equation (41) in the Bloch vector space is given by 

G^{k-u) = G^(k;u) +Y J Go a {hoj)M aP {k;uj)G^{k-uj) (43) 

af3 

The renormalized energies in the mean field approximations are the solutions of the equation 

Y J F^{k)G^{k-u) = l5 ilv (44) 

a 
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Using (44) we find 

Eg{k) = [u- e a (k)]5 au - (1 - 5 au )(V - I)K^ - 

£[(1 - <W)M^ " + (US^ + V(l - SanWwN?]; 

Eg'(fc) = [US a u + V(l - 5 av )]K^- 
N a = iV < ap aa .a pa( j >; 



if. 



a/3 
<T1<T2 



For the degenerate Hubbard model (V = I = J = 0) we get 

EK(k) = [u>- e a (k) - UN?]5 au 



(45) 

(46) 
(47) 

(48) 
(49) 



The spectrum of electronic low-lying excitations without damping follows from the poles of the 
single-particle mean-field GF 



/ En E12 \ / Gqii G012 \ _ / 1 



\ £^2i E22 J V G021 G022 / V 1 



(50) 



Here Go denotes a matrix in the space of band indices. If we put the spin-flip contributions 
equal to zero, i.e. 

E12 = E21 = 
then the quasiparticle spectra are given by 

det|£ii|=0; det\E 22 \ = 

For the multiorbital Hubbard model (5) we find 

GffnH = [u- e a (k) - UN? -VJ2(1- S au )(N? + iVf ) + - M^f]" 1 (51) 

V V 

Finally we turn to the calculation of the damping. To find the damping of the electronic states 
in the general case, one needs to find the matrix elements of the self-energy in (43). Thus we 
have 



G11 G12 
G21 G22 

From this matrix equation we have 



fGon G012V 1 /Mil M X2 \ 

\G02i G022/ V^21 M22) 



-1 



G11 — (G i\ - En) 1 ; G 2 i — (G 02 \ - S 



21 



(52) 



(53) 



G12 = (G 01 2 - S12) 1 ; G22 = (G 02 2 - S22) 1 ; 
where the true self-energy has the form 
En = Mn - E 12 e£M 2 i + \Mi 2 E^ + (M12 - Ei2)E^M 2 2{E22 - M22)- 1 } (M 2 i - #21) (54) 
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The elements of the mass operator matrix M are proportional to the higher-order GF of the 
following form 

{ << a k+paai a p+qvo-2 a qi'V2\ a k+sf3(T3 a rtio-4 a r+stJ,a4 ' J 

For the explicit approximate calculation of the elements of the self-energy it is convenient to 
write down the GFs in (54) in terms of correlation functions by using the well-known spectral 



theorem |53| : 



{ << ak+pacr 1 ap +qV(y2 aq Va2 ^k+sPa^r^a^r+s^ai >> J = 
1 [ +oa duj' r+oo 

7 (exp(/3u/) + 1) / exp(-iu't)dt 



2"7T J-oo <jJ — U)' 



< a 



fc+s/3o- 3 



{t)<l r ^a4,(t) a r+s^a4,(t)\ a k+pa.a 1 &p+qva 2 a q ua ''2 ^ ? ) (^) 



Further insight is gained if we select the suitable relevant "trial" approximation for the correla- 
tion function on the r.h.s. of (55). In this paper we show that the earlier formulations, based on 
the decoupling or/and on diagrammatic methods can be arrived at from our technique but in 
a self- consistent way. Clearly the choice of the relevant trial approximation for the correlation 
function in (55) can be done in a few ways. For example, a reasonable and workable one may 
be the following "pair approximation" ||], which is especially suitable for low density of the 
quasiparticles: 

< a fc+s/3cr(^) a r/i— a(^) a r+Sfi— a{t)Q>k+paoQ>p-\-qv-a a qv-(j > ~ 
< a k+s/3(T(t) a k+pa(T >< a tn-a(t)Q'qv-cr >< a r+S[i-a{t) a p+q U -o > 
+ < a k+sf3o( t ) a 'lv-o >< a rn-v(t) a k+pacr >< «r+s/i-cr (*) a p+g!/-o- > (56) 

Using (56) in (55) we obtain the approximate expression for the self-energy operator in a self- 
consistent form (the self-consistency means that we express approximately the self-energy op- 
erator in terms of the initial GF and, in principle, one can obtain the required solution by a 
suitable iteration procedure): 



pqfiu J 

a/3 i \ . afi / \ vf3 



duJiduj2duJ3 



UJ + 0J\ — L)2 — OJ3 



N(ui,u 2 , w 3 ) [5p+ 3 tT ( w i)% M T ( w 2)Sfc +pTT (w 3 ) + 5fc+ pTT (^ 3 )% T (w2)5p+g TT (wi)] + 

' L ' i ' L ' 2 ' u '' Nfa^usmmm) + (td(it)(tt)] + 



BfBf J 



UJ + UJ\ — U>2 — UJ3 



p j d^d^duis AT(u;i,u;2,W3)[(TI)(IT)(TT) + (TT)(U)(IT)] + 
J UJ + UJ\ - 0J2 — LO3 

e j du^dua Ni^u^sMDaDm) + (U)(U)UT)]) (57) 

J UJ + Ul — Ll>2 — ^3 

where we have used the notations 

N(oji,u 2 ,u 3 ) = [n(uj 2 )n(uj 3 ) + n(u>i)(l - n(oj 2 ) - n(oj 3 ))]; 

gkaa'(u) = --ImGkaa>(u + ie); n(uj) = [exp(/3u) + 

7T 
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Here we present for brevity the explicit expression for a part of the Hamiltonian only without 
the last term. The equations (43) and (57) form a closed self-consistent system of equations 
for the single-electron GF for the Hubbard model, but for the weakly correlated limit only In 
principle, one may use on the r.h.s. of (57) any workable first iteration-step forms of the GFs 
and find a solution by repeated iterations. It is most convenient to choose as the first iteration 
step the following simple one-pole approximation: 

9ka(u) « S(u - e(ka)) (58) 

Then, using (58) in (57), one can get an explicit expression for the self-energy. However, the 
actual explicit calculations will be much more transparent if we confine ourselves to the single- 
orbital Hubbard model in order to discuss more explicitly the reliability of the present approach. 



8 Antiferromagnetic Single-Particle States 

The technique for obtaining of the antiferromagnetic solutions to the correlated fermions on a 
lattice is presented in this section for the single-orbital Hubbard model (10). In general, it can 
be easily applied for the multiorbital extended Hubbard model. 

As discussed above, the self-consistent approach to the calculation of the one-particle properties 
requires the calculation of the folllowing GF 

( <K aiAa h >> << >> ) = G{ir,«) (59) 

\«a i - a \ar T » « ai- a \aj_ a » J 

The equation of motion for the Fourier transform of the GF has the form 

A(im)G(mj; go) = Sij5 af7 i + U « aurTH-alap » (60) 

m 

where 

Him) = ( M ™ " Um) ° (61) 



'i in j 



Using the definition of the irreducible parts (26)the equation of motion can be exactly trans- 
formed to the following form 

J2Mim)G(mj;uj) = 5^5^ + UD ir (ij;co) (62) 

m 

where 

A^im) = ( M ™ - **» " U < >} , s ~ U < ai °°?-° > , ) (63) 
V -U < ai- a al a > [ojdmi -tim-U < n ia >) J 

To calculate the irreducible higher-order GF D tr we have to write the equation of motion for it. 
After introducing the irreducible parts for the operators in the right-hand-side we find 

J2D ir (in;u)A 2 (nj) = C/ 2 £i(u» (64) 

n 

where 

a ,. . ( {ir) « Oiani-ala^nj-a »^ ^ « aiam-alaf^nja \ 

1{IT,UJ) - \M « a^n^n^ »<*■> (*> « a^n^^ »<*> ) (55j 
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Then equation of motion for the GF can be exactly transformed into the following scattering 
equation 

G(ij;u) = G (ij;uj) +^2G (im;uj)P(mn;uj)G (nj;uj) (66) 

mn 

where the generalized mean-field GF Go reads 

^2A 1 (im)G (mj;uj)=5i j (67) 

m 

and the scattering operator P has the form 

p = u2 / W « a ic n^\ajn^ »W W « o^jajy* »W\ 

W ' V (ir) «a l _ (J n iff |a+n i _ (J »M M « a^n^ja+^n^ »W / V ; 

The Dyson equation (41) then will be reduced for the single-band Hubbard model to the following 
form 

G(ij;uj) = Go(ij;uj) + Gp(im; uj)M(mn; u))G(nj; lo) (69) 

mn 

The mass operator M(mn;uj) = U 2 P^ p \mn; oS) describes the inelastic (retarded) part of the 
electron-electron interaction. For the purpose of analogy with the theory of superconductiv- 
ity let us write the Hartree-Fock (elastic) part of the Coulomb mass operator (not included 
in (68)): 

= " Uz^> < <"^ ) > )*- <™> 

To obtain workable expressions for various parts of the mass operator we use the spectral the- 
orem, inverse Fourier transformation and make relevant approximation in the time correlation 
functions. In analogy with the theory of superconductivity the suitable approximation which 
describe the interaction between the charge and spin collective excitations can be written as 









< a ncr( t ) a ma >< n n - a (t)n m _ a 


> 


+ < 4- 




> 


+ < a n - 


-a(t) a m — a ^< a n<r i^) a ri-a if) a ma a m— a 


> 


+ < 


. (t)d m _ cr X Q n _ a (fyQ n — u (£)(l ma Qj m _ cr 


> 


+ < 4- 




> 


+ < a m 




> 



(71) 

The suitable or relevant approximations follow from the concrete physical conditions of the 
problem under consideration. We consider here for illustration the contributions from charge 
and spin collective degrees of freedom. We get 

n U 2 f+°° , , ctg^+tg^ 
M(i 3] lu) = — / duj 1( luJ2- 2 



(( 



2VT 2 7_oo LU — UJl — LU2 

f Im « ni_ a \rij^ a » Wl Im « a i(T \a^ a >> W2 Im « ni- a \n ja » Ul Im « a ia \a^_ (T >>^ 2 

Im « n ia \nj- a » Ul Im « ai- a \aj a >>a, 2 Im « n ia \n j(J >> Wl Im « ai_ CT |a+_ CT >>a, 2 

flm« S^ a \Sj » u>1 Im« ai- a \a+_ a » W2 Im « S^ a \Sj a » u>1 Im « ai- a \aj a » u>2 

V Im « S°\SJ » W1 Im « a ia \a+_ a » W2 Im « S°\Sy a » Wl Im « a ia \a+ a » U)2 
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It shows that it is possible to do all calculations in the localized Wannier basis as we did in 
deriving the equations for the strong coupling superconductivity in transition metals [19]. This 
has great advantage for consideration of disordered transition metal alloys. 
As for the translationally invariant crystal with broken symmetry the following special Fourier 
transform should be performed for the generalized mean-field GF Go(ij;to) (67) 

G$(ij;u) =J2 ex P[ik(Ri-R j )]G 1 1 (k;u) (73) 

k 

Gl 2 {ij;to) = J2 ex P [ikRi - i{k + Q)R j ]G 1 2 {k;uj) 

k 

Gf{ij;u) = J2 ex P [K k + Q) R i ~ ikRj]Gf(k;u) 

k 

Gf{ij-uj) = ]Texp [i(k + Q)(Ri - R^Gfik;^) 
k 

The result of this transformation is then 



u -Ef F (k + Q) A n (k) 



G °-{Gf Gf)~ (oj-Er(k))(co-E^(k)) (74) 
where 

E? F = e(k) + U <n a > (75) 
A a - a {k) = U ^2iexp(ikRi) < a^af^ > 

W - ( E " F{k) + T k + Q) * i C fF(k) ~ T { Q) r ) 

It is evident that one can define the Overhauser's angle 6k 

cos 2 « A n (fc)A iT (fc) 

k ~ {u-E?r{k)F + ± n (k)A u (k) [7b) 

In Overhauser's notations A^(k) = A^(k) = A . The self-consistent set of equations for 
determining of the SDW (or "gap") order parameter A , chemical potential [i and averaged 
moment < s z > is 



k 

<s z > ^ 



A = [//iV^<a+ +Qi a fcT > (77) 

k 
k 

n = N-^^iEf^ik) + n(Ei IF (k)) 
k 

The above expressions were derived for correlated itinerant fermions on a lattice within the 
Hubbard model and for finite temperatures. These equations were also deduced in previous 
papers in the course of their analysis. Here we deduced it by using more sophisticated arguments 
of the IGFs method in complete analogy with our description of the Heisenberg antiferromagnet 
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at finite temperatures However, the self-consistent system of equations (69), (72) for 

determining the quasiparticle spectra with damping is not as obvious as generalization as the 
equations (77). This is intrinsically the many-body manifestation of the correlation effects of 
itinerant fermions on a lattice and shows clearly the advantage of the present approach. 
To confirm this, the explicit calculation of the damping should be performed. The natural way to 
tackle this program would then to look at the calculations of the collective GFs of the generalized 
spin (and charge) susceptibilities in (72) but this deserves of separate consideration. Again this 
problem bears close similarity to the paramagnetic Hubbard model and the antiferromagnetic 
Heisenberg model and it can be argued that this effect of interference of single-particle and 
collective modes of excitations should be considered carefully. 



9 Discussion 



We have been concerned in this paper with establishing the essence of single-particle excitations 
of correlated lattice fermions, rather than with their detailed properties. We have considered 
the single- and multiband Hubbard model but the calculational details were mainly presented 
for the single-band Hubbard model where the appropriate concepts are easier to demonstrate. 
We have considered a general family of symmetry broken solutions for itinerant lattice fermions, 
identifying the type of ordered states and then derived explicitly the functional of generalized 
mean fields and the self-consistent set of equations which describe the quasiparticle spectra and 
their damping in the most general way. While such generality is not so obvious in all applica- 
tions, it is highly desirable in treatments of such complicated problems as the competition and 
interplay of antiferromagnetism and superconductivity, heavy fermions and antiferromagnetism 
etc., because of the non-trivial character of coupled equations which occur there. Both these 
problems are subject of current but independent research. 

Another development of the present approach is the consideration of the itinerant antiferromag- 
netism of highly correlated lattice fermions when U is very large but finite. Like the weakly- 
coupled case described in this paper, the symmetry broken approach will work, but matters 
are complicated by the necessity of constructing the more extended algebra of relevant opera- 
tors H . This idea has been carried out for the paramagnetic solution of the single-band Hubbard 
model [ 10 1 . It would be interesting to understand on a deeper level the relationship between 
Mott-Hubbard metal-insulator transition and various ordered magnetic states within the Hub- 
bard model. 

In conclusion, we have demonstrated that the Irreducible Green's Functions approach is a work- 
able and efficient scheme for the consistent description of the correlated fermions on a lattice 
at finite temperatures and that it can be generalized naturally to include the symmetry broken 
concept. 
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